• 



COMPUTER SYSTEM 
100 



MEMORY 112 




COMPUTATIONAL 






CODE FOR 






INTERVALS 






114 




I 




INTERVAL VALUES 






115 







fn 

ft==s 








BRII 
- , 1( 


DGE 
)6 




FIG. 1 



PERIPHERAL BUS 
110 




COMPILER FOR 
INTERVAL CODE 
204 




ARITHMETIC 
UNIT FOR 
INTERVAL 
COMPUTATIONS 
104 



FIG. 2 



* 



• 




a 



-4 

til 











ING 


ER3 






0 
_j 


IUMI 


LL 


z 


H 


H 


CO 




a: 


O 


LL 


CL 



CO 

0 





DATING 
BER 306 




FLOATING 
NUMBER 314 




DATING 
BER 316 


NID FL 
NUM 






NID FL 
NUM 


0 i 








0 i 


HI O 




%. O 
ll. a. 




ill O 


a> Q. 










Q co 
Z O 
LU ^" 



a 

H 
LU 



CO 
3 



So 



a: 
o 

LL 
01 
LU 
CL 



Ll 
O 

O o 

LU * 

CO w 

LU t= 

(T ? 
CL 
LU 



6 



o 
< 

LU < 

CO > 

LU Ll 

* £ g 
CL t= tJ- 
LU ? 

* LL 

w O 

CD 
o 

LU 



H 

0£ o 
< 2 

CO 



□ 



IK 

VI 

x 

VI 
X 



IK 

XI 



VI 
VI 
■3* 



I?S 
>%l 

1 ■ 



"I I" 



i;\ 
+ 



>>i 

i 



IK 


ik 


<- 


<- 




1^ 


+ 


1 


XI 


XI 


-> 

1 1 


-> 


II 


II 


+ 


1 


X 


X 



X 
IK 

>^l 
X 
IK 

X 

XI 

>>\ 

X 

XI 

3 



1^ 

X 
IK 

X 
IK 

rs 

X 
XI 

>>\ 

X 

XI 

^5" 



ii 

X 

X 



O 



IK 

IK 
1^ 

"Si 

"Si 
<- 

B 

IK 
>>\ 

IK 
1^ 

"Si 
>,l 

XI 



X 



w 
o 

Ul 
X 



in 
CD 




C START J 



RECEIVE REPRESENTATION OF 
NONLINEAR EQUATION /x)=0 
602 



RECEIVE REPRESENTATION OF 
INITIAL INTERVAL X 
604 



SYMBOLICALLY MANIPULATE 
7(x)=0 TO SOLVE FOR THE TERM 
g(x')=/i(x), WHEREIN g(x') CAN BE 
ANALYTICALLY INVERTED 
606 



SUBSTITUTE X INTO g(X')=h(X) 
608 



SOLVE FORX'=g-'[/z(X)] 
610 





INTERSECT X' WITH X TO 
PRODUCE NEW INTERVAL X + 
612 
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PERFORM INTERVAL NEWTON 
STEP ON/x) AND X TO PRODUCE 
NEW INTERVAL X + 
616 
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EVALUATE/AT CENTER OF EACH 
INITIAL BOX IN L, AND USE 
RESULT TO \1?D ATE f bar 
701 



FOR EACH INITIAL BOX IN L 
EVALUATE/X). DELETE ANY 
BOX FOR WHICH THE LOWER 
BOUND OF XX) IS >f_bar 
702 



I 



IF L, IS EMPTY, GO TO STEP 738. 
OTHERWISE, FIND THE BOX X IN 
L, FOR WHICH THE LOWER 
BOUND OF/X) IS SMALLEST. 
DENOTE THIS BOX AS X(». 
DELETE XW FROM L r 
703 



EVALUATE/AT CENTER OF X 
AND USE RESULT TO UPDATE 
fjbar 

(SKIP THIS STEP IF ALREADY 
DONE FOR BOX IN STEP 701). 
704 



IF TERM CONSISTENCY HAS BEEN 
APPLIED n TIMES IN STEP 706 
WITHOUT HAVING DONE STEP 
710, GO TO STEP 709. 
705 



APPLY TERM CONSISTENCY TO 
Ax) zf_bar. IF RESULT IS EMPTY, 
DELETE X AND GO TO STEP 703. 
706 



IF w(X) < e x AND w(/[X)) < e p PUT X 
IN LIST L 2 AND GO TO STEP 703. 
707 



IF BOX WAS SUFFICIENTLY 
REDUCED, PUT RESULT IN L, AND 
GO TO STEP 703. 
708 



IF TERM CONSISTENCY HAS BEEN 

APPLIED n TIMES TO 
COMPONENTS OF THE GRADIENT 
WITHOUT HAVING APPLIED A 
NEWTON STEP, GO TO STEP 722 
709 



APPLY TERM CONSISTENCY TO 
g,.(x)=0 (7=1, n) FOR EACH 
COMPONENT g ( (x) OF THE 
GRADIENT OF/fr)-' IN DOING SO, 
BOUND g t OVER THE RESULTING 

BOX. IF RESULT IS EMPTY, 
DELETE C AND GO TO STEP 703. 
710 



USE THE CENTER OF THE 
BOUNDS OBTAINED FROM 
WORKING ON g ( FROM STEP 710 
TO DO A LINE SEARCH AND 
UPDATE/J>ar. 
711 
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USING BOUNDS FROM WORKING 
ON g.(X) (/=1, ri) OBTAINED IN 
STEP 710, APPLY THE LINEAR 
METHOD TO THE RELATION 
Ax)sf_bar. 
712 



REPEAT STEP 707 
713 



IF CURRENT BOX X IS A 
SUFFICIENTLY REDUCED VERSION 
OF XV\ PUT X IN L l AND GO TO 
STEP 703. 
714 



IFjV„=4, GO TO STEP 720. 
715 



APPLY TERM CONSISTENCY TO 
//..(X) ;> 0 (z'=l, n). IF RESULT IS 
EMPTY, GO TO STEP 703. 
716 



IF BOX CHANGED IN STEP 716, SET 
N H =0 AND GO TO STEP 718. 
OTHERWISE, IF BOX IS 
UNCHANGED IN STEP 716, REPLACE 
N H BY N H +\ AND GO TO STEP 720. 
717 



REPEAT STEP 707. 
718 



REPEAT STEP 714. 
719 



REPEAT STEPS 706 THROUGH 719 

USING BOX CONSISTENCY IN 
PLACE OF TERM CONSISTENCY. 

HOWEVER, SKIP STEPS 71 1 
THROUGH 714. IN STEP 710, OMIT 
THE PROCESS OF OBTAINING 
BOUNDS ON g ( .(X). 
720 



IF w(X) > (w, + w R )/2, 
GO TO STEP 732. 
721 



COMPUTE J(x,X) OF THE GRADIENT 
g. COMPUTE AND APPROXIMATE 
INVERSE B OF THE CENTER OF 
J(x,X) AND THE MATRIX 
M(x,X)=BJ(x,X). 
722 



IF M(x,X) IS REGULAR, FIND THE 
HULL OF THE SOLUTION SET OF 
THE LINEAR SYSTEM DETERMINED 

IN STEP 722. IFM(x,X)IS 
IRREGULAR, APPLY ONE PASS OF 
GAUSS-SEIDEL METHOD TO THE 
LINEAR SYSTEM. UPDATE w, OR w R . 
IF RESULT OF THE NEWTON STEP IS 

EMPTY, GO TO STEP 703. IF 
INTERVAL NEWTON STEP PROVES 
EXISTENCE OF A SOLUTION IN X, 
RECORD THIS INFORMATION. 
723 
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USE GRADIENT VALUE g(x) AND B TO 
COMPUTE THE POINT g=x-Bg(x). USE 
THE VALUE 0¥f(y) TO UPDATE/J>ar. 
724 



USE QUADRATIC METHOD TO SOLVE 
/x) <f_bar. 
725 



REPEAT STEP 707. 
726 



REPEAT STEP 714. 
727 



USING B, ANALYTICALLY DETERMINE 
Bg(x). APPLY TERM CONSISTENCY TO 
SOLVE THE / th COMPONENT OF Bg(x)=0 
FOR THE i* VARIABLE x j (i=\, n). IF 
THIS PROCEDURE PROVES THE 
EXISTENCE OF A SOLUTION IN X, 
RECORD THE INFORMATION. 
728 



REPEAT STEP 707. 
729 



REPEAT STEP 714! 
730 
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APPLY BOX CONSISTENCY TO THE / th COMPONENT OF 
Bg(x) FOR THE z th VARIABLE FOR i=\, n. 
731 



REPEAT STEP 707. 
732 



IF WIDTH OF BOX WAS REDUCED BY A FACTOR OF AT 
LEAST 8 IN THE NEWTON STEP, GO TO STEP 722. 

733 



REPEAT STEP 714. 
734 



MERGE ANY OVERLAPPING GAPS IN COMPONENTS OF 
X. IF ANY GAPS ARE SUITABLE FOR SPLITTING X, GO 
TO STEP 733. OTHERWISE, GO TO STEP 737. 

735 



SELECT THREE GAPS AND SPLIT X INTO SUB BOXES 
THEN GO TO STEP 703. 
736 



SPLIT THE THREE COMPONENTS OF X FOR WHICH D. 
(/=1, n) IS LARGEST. SPLIT EACH SUCH COMPONENT, 
THEN GO TO STEP 703. 
737 



DELETE ANY BOX X FROM LIST L 2 FOR WHICH THE 
LOWER BOUND OF/X) IS LESS THAN/J>ar. DENOTE 
THE REMAINING BOXES X<'\ X®. DETERMINE LOWER 
BOUND FOR GLOBAL MINIMUM. 
738 



TERMINATE. 
739 
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